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In this paper, there are two contributions: The first contribution is to design a 
robust cascade P-PI controller to control the speed and position of the 
permanent magnet DC motor (PMDC). The second contribution is to use 
three methods to tuning the parameter values for this cascade controller by 
making a comparison between them to obtain the best results to ensure 
accurate tracking trajectory on the axis to reach the desired position. These 
methods are the classical method (CM) and it requires some assumptions, 
the genetic algorithm (GA), and the particle swarm optimization algorithm 
(PSO). The simulation results show the system becomes unstable after 
applying the load when using the classical method because it assumes 
cancellation of the load effect. Also, an overshoot of about 3.763% is 
observed, and a deviation from the desired position of about 12.03 degrees is 
observed when using the GA algorithm, while no deviation or overshoot is 
observed when using the PSO algorithm. Therefore, the PSO algorithm has 


superiority as compared to the other two methods in improving the 
performance of the PMDC motor by extracting the best parameters for the 
cascade P-PI controller to reach the desired position at a regular speed. 
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1. INTRODUCTION 

A permanent magnet DC motor (PMDC) is a simple type of DC motor where field windings have 
been replaced by permanent magnets. This motor works with the same basic principle as a shunt-connected 
motor but the difference between them is that the permanent magnet generates the required flux instead of the 
field windings [1], [2]. The power of these motors is less than that of the shunt-connected motors because the 
flux generated by the permanent magnet is less than the flux generated by the field windings. Therefore, most 
types of small DC motors are PMDC motors and these motors operate at high speed and low torque [1], [3]. 
There are many applications for the PMDC motor, such as moving windows in cars and front areas of cars. It 
is also used in children's toys and household appliances like food mixers and others [4], [5] and its important 
use in computer numerical control (CNC) machines as an actuator, electric vehicles, and robotics [6], [7]. 

The main problem statement in this paper is how to get the best parameter values for the 
proportional integral derivative (PID) cascade controller to give the most accurate results for tracking the 
trajectory of the reference position to reaching the desired position at a regular speed. To address this 
problem, a comparison was made between the three methods classical method (CM), genetic algorithm (GA), 
and particle swarm optimization algorithm (PSO) to extract the best parameters of the cascading controller, in 
addition to using the objective function, integral time absolute error (ITAE) to reduce the error between the 
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input and output of this system. The cascade P-PI controller used in this paper consists of three controllers: 
the current controller as an inner loop, the speed controller, and the position controller as outer loops [8]. P, 
PI, and PI controllers are used for position, speed, and current respectively [9]. 

The purpose of using the cascade P-PI controller is due to many reasons, the most important of 
which is robustness to reduce or reject external disturbances and back to the steady-state in a short time [10]- 
[12]. There are several methods to tuning the parameters of the cascade P-PI controllers such as the ziegler- 
nichols (Z-N) method, cohen-coon method, neural network, and fuzzy logic [13], [14]. Particle swarm 
optimization (PSO) algorithm [15], [16] and genetic algorithm [17]. 

To compare this work with other researcher's work in the same field, the following literature review 
was made: In 2014, Mustafa et al. [18] presented a brushless DC (BLDC) motor speed control system using 
GA. In 2015, Taha et al. [9] used three methods to control the cascade control system. In 2018, Wisam et al. 
[19] presented a system for controlling PMDC speed using GA and direct search (DS) algorithms. In 2019, 
Fadhel et al. [20] used a fractional PID controller to control PMDC speed based on PSO. In 2021, Ahmed et 
al. [21] presented a system to control the position and speed of a servo motor. This paper is organized as: The 
second section contains the mathematical model of the PMDC motor and explaining the general structure of 
the system, in the third section, the three tuning methods are explained with the objective function (ITAE). 
The fourth section contains the results and comparison, and the final section contains the conclusion. 


2. MATHEMATICAL MODEL AND GENERAL STRUCTURE 
2.1. Mathematical model of PMDC motor 

Figure 1 shows the equivalent circuit of a PMDC motor consisting of an armature resistance (Ra) 
and inductance (L,) connected in series. The back emf (E,) is generated when the flux lines generated by the 
permanent magnet are cutting and its direction opposite to the direction of the applied voltage. While the 
mechanical part consists of coefficient of friction ( Bm) and moment of inertia (Jm). In addition to other 
parameters are the back emf constant (K,) and torque constant (K,). All components or parameters of this 
PMDC motor and their values are shown in Table | and Figure 2 shows the block diagram of PMDC motor. 


Figure 1. Equivalent circuit of PMDC motor [22] 


Table 1. PMDC parameters [9] 


Motor parameters Value 
Torque constant K,=2.35 Nm/A 
Armature inductance La=2.61*10'°H 
Armature resistance Ra=2.61 Q 


Inertia of the motor Jm=0.068 kg.m? 
Friction constant Bm=0.008 Nms/rad 
Back emf constant Kv=2.35 Vs/rad 
Nominal load T,=17.6 Nm 
Nominal voltage Va=230 v 


Figure 2. The block diagram of PMDC motor [22] 
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The electrical and the mechanical equations corresponding to block diagram from Figure 2 are [22]. 


Va (t) = ea (t) + Ra ia (Ò + La Sia (t) (1) 
ea (t) = ke Om (t) (2) 
Tm (t) — T, = Sm ge Om (t) + Bm (t) (3) 
Tm (t) = Keia (Ò (4) 


By using Laplace transformation for the (1)-(4) we obtain; 


Va (Ss) = Eq(s) + Ra la(s) + SLgla(s) (5) 
Ea(s) = ke @m(S) (6) 
Tm (s) =T, = SJm@m(S) + Bn@m(S) (7) 
Tm (S) = K; la(s) (8) 


The overall transfer functions are defined for speed and position control of PMDC motor, respectively. 


wms) = Ke (9) 
Va (S) J LqS2 +(JRq +BLg )S+BRq + ke” 


a(s) _ Ke 
Va (S) J LaS? +(JRa +BLq)S? +( BRq +ke")S 


(10) 


Where ke is equal to K, [23] and 


Ra = armature resistance (Q) 

La = armayure inductance (H) 

Ea = electromotive force or back emf (v) 
Q(s) = actual position (rad) 

Va = nominal voltage (v) 

Jin = moment inertia (kg.m”) 

TL = nominal load torque (Nm) 

Kı = torque constant (Nm/A) 

Ky = back emf constant (v.sec/rad) 
Bn = friction coefficient (Nm.s/rad) 
Om = motor velocity (rad/sec) 


2.2. General structure of the system 

The general structure of the system consists of a cascade P-PI controller as shown in Figure 3. This 
controller consists of three (P, PI) controllers for current, speed, and position. The output of the position 
controller represents the reference speed, the output of the speed controller represents the reference current, 
and the output of the current controller represents the control voltage (VC). 
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L P ref BU oe 


ITAE 


Figure 3. General structure of cascade PID control system for PMDC 
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An advantage of the cascade control system, in addition to the other advantages, is the ability to set 
limits in order to protect the PMDC motor and the power electronic converter. In this paper, limits are placed 
on the speed reference by an amount not exceeding the motor's rated speed. As well as the limit, placed on 
the reference voltage exiting the PI current not exceeding the motor's applied voltage [23]. 


3. TUNING METHODS 
3.1. Classical method (CM) 

Figure 4 represents the first control loop or the inner control loop. Since this method assumes some 
assumptions to simplify the circuit, the effect of the torque load is neglected, and the effect of E, is neglected 
because the value of J,, is high [9], [23]. As a result, the simplified current control loop as shown in Figure 5. 


T,(S) 


(S) He 


Ea (S) 
< 


Figure 5. Simplified inner current control loop 


From Figure 5 the transfer function of the current control loop, G; 9;(s), and it can be written like in 
(11) where, ki; is the integral gain, k,, is the proportional gain of the controller, and the electrical time 
constant Te, is calculated using (12). 


1 
_ ki S / Ra 
Gi oS) == (1+ Ku fe -). ren) (1) 
p 
lig 
te = ot (12) 


The (13) is used to cancel the motor pole, which is formulated in (11). 


Ku = + (13) 


kp Te 


The cancellation of the pole in the motor transfer function can be illustrated as: 


1+5STe HESE 


1 1 
Gror = a (iz ): ( Hes.) OR Go = U (1 + sty) (42) 
Kp 


Finally, in the open-loop transfer function, the bandwidth (crossover) frequency (wci) It is represented by 
(14). 
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Also, the crossover frequency of the current open-loop can be calculated from the following 
relationship wci = 2mfci and its value about ten times smaller than the switching frequency for DC-DC 
converter [9]. So, the parameter values of the PI current controller can be found from (13), (14). To find the 
speed loop parameters it is assumed that the closed current loop is ideal for design purposes and is 
represented by the unity [23] as shown in Figure 6. 


Figure 6. Outer speed loop 


From Figure 6, the transfer function of this loop Gs ,;(s), can be represented in (15). 


Grols) = 2 (14 5 |» (2m 15 
soils) = -> Hs] Cree, (15) 
p 
Where 
Tm= mechanical time constant = da (16) 


Just like the inner current control loop the pole in the mechanical part of the motor will be canceled as show 
in (17). Here, kis is the integral gain and, kps is the proportional gain of the controller. 
kis 1 


=— (17) 


kps Tm 


The bandwidth (crossover) frequency of the speed control loop, wcs can be chosen to be ten times lower than 
oci and it can be represented by (18). 


— kis-kt 
ocs = -57 (18) 


Thus, kiş and k„s can be calculated using (16) and (18). 

Finally, the last outer loop is the position control loop it can be illustrated in Figure 7. To find the 
k_pP parameter of the position loop it is assumed that the speed loop is perfect and is represented by unity 
[23] as shown in Figure 7 the open-loop transfer function for position control is given in the (19). The value 
of the position parameter k_pP can be calculated from (20) where the bandwidth frequency (wps) is chosen 
to be ten times smaller than (acs). 


Figure 7. Outer position loop 
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Gp oi(s) = E (19) 


P 
S 
æps = kpp (20) 


3.2. Particle swarm optimization (PSO) 

Particle swarm algorithm (PSO) is proposed by Kennedy and Eberhart and it was modified to improve 
its performance by adding a new parameter called inertia weight [24]. The PSO algorithm bases on swarm 
intelligence techniques for the observation of social behavior of moving organisms such as a gathering of fish or 
birds. This algorithm is related to the computational method to solve the problem, for example, bird flocks aim 
to find eating behaviors and use repeated steps to reach the best solutions. This means that the candidate 
solutions are particles that move in the search space based on a specific formula above the particle position. Each 
particle's movement is affected by its local value, and its objective is to reach the best-known positions in the 
search space by updating better positions found by other particles. 

The algorithm is begun by establishing the starting position and speed vectors. At each iteration, the 
best value is determined by evaluating position and speed vectors. Every particle has variables and dimensions, 
and these variables are problems that need to be solved. If the problem consists of five different variables, the 
particles’ dimension should be chosen as five. Each particle’s best value is called a local best value and recorded 
into the P best matrix. After each iteration, the best value for each particle is updated if the best new value is 
found to control each current particle and the previous positions. Besides, the position is affirmed as the global 
best after controlling for the best matrix values at each iteration which continues until it reaches the specified 
number, as the last updated value for the best position represents the optimal value. This algorithm depends 
mainly on finding the position of each particle with the best local value, as well as finding the best general 
swarm position in each iteration. The position and speed are updated at each iteration based on (21) and (22) 
[20], [25]. 


Vij (t + 1) = W. Vij (t) + rıcı [Pbest; ;Ct) =, Xij (t)| 


+ T2C2 [Gbest; j (t) = Xij (t)] (21) 
Xit +1) = Vij +1 + Xy) (22) 

Where 

i = Particle index 

j = Dimension number 

t = Iteration 

Vi; (t) = Current velocity 

Xij (t) = Current position 


V,;(t + 1) = Velocity updated or new velocity 
Xi j(t + 1) = Position updated or new position 


W = Weighted inertia its value between 0 and 1 
Ti’ = Random coefficients there values between 0 and 1 
C1 C2 = Acceleration coefficients there values between 0 and 2 


In Table 2 are the parameters of the PSO used in this work and Figure 8 shows the flowchart of the 
particle swarm optimization algorithm. 


Table 2. PSO parameters 


PSO Parameters Value 
Iteration 100 
Swarm size 20 
No. dimension 5 
Weighted inertia 0.9 
Cl 2 
C2 1.5 
LU, UB 0, 300 


3.3. Genetic algorithm 
In 1970, Holland introduced the genetic algorithm (GA). This method is based on the principle of 
evolution and natural selection. As an optimizer, this algorithm begins by selecting a specific number of 
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chromosomes from a primary set, and each chromosome represents a solution to a specific problem and its 
performance is evaluated by the fitness function [18]. The present parents produce new offspring called the 
new generation, as this new generation is supposed to reach a better solution than previous parents, due to the 
process of evolution. This mechanism is repeated for several successive generations until one of them obtains 
the optimal solution for the system. Since the process is associated with many phases, Figure 9 shows a 


flowchart of the genetic algorithm [26]. Table 3 represents the parameters of the GA that were used in this 
work. 


Initia lized 
Calculate the 
objective function 
Determine Pbest 
Determine Gbest 
Update velocity 
and position 


Check the stop 
condition 


CREATE 
POPULATION 
FIRST 
GENERATION 
CALCULATE 
FITNESS VALUE 
SELECTION 


CROSSOVER 
MUTATION 


Increase 
Iteration number 


OR MIN. P IN, 
REACHED 


Figure 8. PSO flowchart [27] Figure 9. Genetic algorithm process flowchart [26] 


Table 3. GA parameters 


GA Parameters Value 
Generation 100 
Population size 20 
No. dimension 5 
Crossover 0.6 
LU, UB 0, 300 


When designing any control system for a particular physical process, certain criteria must be chosen 
to provide the best performance. These criteria are known as objective functions. These functions reduce the 
error between the input and output and have several types. In this work, a performance indicator was used to 
reduce the error in the position signal called an integral absolute time error (ITAE) and can be represented by 
(23). 


ITAE = J, t|e(t)|dt (23) 


4. RESULTS AND COMPARISON 

The simulation of the system was performed using MATLAB 2019, and the system was tested after 
entering the cascade P-PI parameters extracted from the three tuning methods. Table 4 shows the values of 
these parameters extracted from each method. The system has also been tested in other cases, such as the load 
condition, no-load, and single or multiple reference position to ensure the PMDC motor rotates 360 degrees 
or several revolutions. Table 5 shows the performance parameters for the system at load case, which were 
obtained from each tuning method. 
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Table 4. PID parameters values 


PID Parameters CM GA PSO 
KP. Position 125.6637 244.3101 173.1189 
KP. Speed 36.3623 98.7314 269.2364 

KI. Speed 4.2779 291.317 3.6981 
KP. current 32.7982 101.5191 216.5567 
KI. current 32798 268.9725 23.8736 


Table 5. The values of the performance criteria for each tuning method 


Performance criteria CM GA PSO 
Rise time 2.129 sec 2.033 sec 2.053 sec 
Settling time 0.3 sec 0.2 sec 0.2 sec 
overshoot 7.558% 3.763% 0 


Through the simulation results attached in the appendix, and for the purpose of obtaining the best 
results, the following comparison was made between the three tuning methods: 

a) At no-load case, high overshoot is observed in speed when using the CM method as shown in 
Figures 10(a) and 10(b), and after applying the full load at the fifth second, the system becomes unstable 
because this method assumes neglecting the load effect to simplify the system as shown in Figures 1 1(a) 
and 11(b). 

b) When using the GA algorithm and at the no-load case, an overshoot was observed about 3.646 %, and a 
deviation of 12.03° was observed from the desired position, as well as a high overshoot was observed in 
the speed as shown in Figures 12(a) and 12(b). After applying full load, a very small deviation of about 
0.0003 rad was observed. This is an indication of the robustness of the cascade P-PI controller in 
rejecting external disturbances as shown in Figures 13(a) and 13(b). 


40 r r r r 
Ref 


overshoot in speed 
at no load case 


Position [rad] 
Speed [rad/sec] 


Time [sec] Time [sec] 


(a) (b) 


Figure 10. These figures are: (a) position control at no load case, and (b) speed control at no load case 
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Figure 11. These figures are: (a) position control at load case, and (b) speed control at load case 
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Figure 12. These figures are: (a) position control at no load case and (b) speed control at no load case 
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Figure 13. These figures are: (a) position control at load case, and (b) speed control at load case 


When using the PSO algorithm, there was no deviation from the position or overshoot in speed or 
position in both cases load and no-load in addition to the accurate tracking trajectory to reach the 
position with high accuracy as shown in Figures 14(a), (b) and Figures 15(a), (b). The negative speed 
means that the PMDC motor is rotating in the opposite direction. Therefore, the simulation results and 


the comparison proved that the particle swarm optimization algorithm gave better results than the other 
two methods. 


7 1 id | | | Ref 
—— PSO 
6f Ref |] FSO 
5 — 20} 
3 
4 £2 
3 
3} £ 
3 
2+ | pr] 
2 
N 
1 -20 
0 
-1 -40 
0 2 4 6 8 10 0 2 4 6 8 10 
Time [sec] Time [sec] 
(a) (b) 


Figure 14. These figures are: (a) position control at no load case, and (b) speed control at no load case 
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Figure 15. These figures are: (a) position control at load case, and (b) speed control at load case 
CONCLUSION 


This paper presents a comparison between two intelligence methods which are PSO and GA 


algorithms in addition to the traditional method CM. These methods are used to extract the best parameters of 


the 


cascade P-PI controller to control the speed and position of the PMDC motor to reach the desired 


position. The simulation results prove the effectiveness of the PSO algorithm in terms of no overshoot as 
compared with CM and GA methods which have an overshoot about 7.558% and 3.763%, respectively. Also, 
there is no deviation in position when using the PSO algorithm while a deviation about 12.03 degrees is 
observed when using GA algorithm. So, the PSO algorithm has superiority in improving the overall 
performance of the PMDC motor by extracting the best parameters for the cascade P-PI controller to reach 
the desired position as compared to the other two methods. 
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